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We investigated the charge distribution in Pri+a;Ba2-a;Cu306+H-crystals at hquid-He temperature 
by means of ^'^'^^Cu- and ^''^Pr-NMR and NQR. The electric field gradients determined for the 
different oxygen coordinations of Cu(l) on the chains confirm the model that the hole states are 
localized in the 02p7r-orbitals of the Cu02-planes. The intensity and line-shape analysis of the 
Cu(l)- and Pr-resonances in the oxygen doping series (a; ~ 0, 1/ = . . . 1) and in the Pr/Ba solid- 
solution series (at y = 1) allows us to assign the Pr-resonance to Pr at RE-sites. Therefore, we 
can investigate the charge distribution in the Cu02-Pr-Cu02-layers by crystal field analysis of the 
Pr-signal. To consistently describe our results with neutron scattering data we propose that two 
Pr-states are present in the RE-layer, and ascribe them to Pr with and one without a hole localized 
in the oxygen coordination shell. NMR detects only the former, with a virtually nonmagnetic singlet 
ground state, while space-integral techniques are dominated by the latter, due to a quasi-doublet 
ground state and antiferromagnetism of its Pr-moments at low temperature. 
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I. INTRODUCTION 

PrBa2Cu307 is the one insulating member of 
the isostructur|al series REBa2Cu307 of high-Tc 
superconductorsEJ (RE=Rare Earth, Y, La, abbreviated 
below RE 123). Simultaneously the magnetic properties 
of Pr differ dramatically from the reasonably well un- 
derstood 4/— magnetism of all other rare earths in this 
structure. This extraordinary behaviour of Pr attracted 
considerable interest from experimental and theoretical 
physics throughout the last decade because it might shed 
some light on the microscopic origin of high-Tc supercon- 
ductivity in the cuprates. From the point of view of ap- 
plied physics it is important to understand the physical 
properties of the material because of its potential use in 
superconductor-insulator heterostructures. 

There is a general consensus that the structural key 
element of the cuprate superconductors are the CUO2- 
planes. The undoped Cu02-plane is an antiferromagnetic 
charge transfer insulator: Photoemission spectroscopy 
shows that the on-site correlation energy of the Cu- 
dx2-y2 electrons is large enough to push the lower Hub- 
bard band below the 2p-oxygen bandH. Upon hole doping 
an insulator-metal transition (MIT) takes place, with a 
superconducting ground state in a certain range of hole 
concentrations. The dependence of the superconduct- 
ing transition temperature Tc on the hole concentration 
rifi in the Cu02-plane follows a universal law close to 
Tc/Tcmax = hn (("'i ~ nopt)/^nY, where this has posi- 
tive solutionajQ (note that other dependencies have been 
proposed, e.g. the trapezoid shape in refE). In con- 



trast to the large differences in the maximum transition 
temperature Tcmax between cuprates, the variations of 
the optimum hole concentration riopt ~ 0.3/CuO2-unit 
and the threshold concentration A„ for superconductiv- 
ity appear to be rather small. The fact that Prl23 is an 
antiferromagnetic insulator instead of metallic and super- 
conducting like the other members of the RE123 series 
may then be due to either a suppression of the doping, 
or a detail of the electronic structure lowering Tcmax dra- 
matically. 

There is no obvious and generally accepted parameter 
connected with Tcmax of the cuprates. In RE123 one finds 
nearly independent of the RE Tcmax = 93±4 K. Magnetic 
pair breaking by the Pr-4/ moments, enhanced by a hy- 
bridization with the bands in the Cu02-planes, has been 
invoked as a mechanism to suppress superconductivity 
by lowering Tcmax- A phenomenological combination of 
Abrikosov-Gorkov pair breaking and hole depletion can 
successfully describe the dependence of the supero 
ducting Tc on the Pr-concentration in Pr2REi_^123J 
and it has been argued that the exceptionally high or- 
dering temperature of the Pr-sublattice, Tjy — 12 .pj. 18 K 
compared to 2.1 K for Gd, also supports this viewE How- 
ever, the results cannot be described by pair breaking 
alone, it is necessary to take hole depletion into account. 
Otherwise Pr would not suppress superconductivity and 
the MIT at the same critical Pr-concentration, an 
would drop linearly at small x, which is not observed 
Furthermore, there is evidence against pair breaking,ll2 
and it is possible to model Tc{z) without magnetic pair 
breaking under the assumption that the influence of Pr 
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on the doping depends on its concentrationJatil Even if 
one agrees that pair breaking is indicated experimentally 
at least near the critical Pr-concentration, it is not ob- 
vious that a Pr-moment is involved since near the MIT 
localized Cu-nioments in the Cu02-layers may be candi- 
dates for pair breakers as well. 

Therefore, the current attempts to explain the role of 
Pr focus on the suppression of the hole doping. There 
is clear experimental evidence that Pr does not sim- 
ply fill the-hpie states with its 4/-electrons. Various 
technkupsJiSiij including inelastic neutron .saattcring 
(INS)lla't3 and our own previous NMR work,tZHl3 show 
that the valence of Pr is close to three, the same as 
for the other RE,|-a|nd stays constant independent of the 
oxygen content xilB Furthermore, spectroscopic methods 
sensitive to the density of 0-2p hole states-detect similar 
amounts in Prl23 as in the other RE123.EIeI 

The by now generally accepted model for the suppres- 
sion of hole doping was proposed by Fehrenbacher and 
Rice. In their model the holes become localized in the 
antibonding 0-2p7r-orbitals of the planes.E3 These states 
are favourable to the 2pfT-states populated in the absence 
of Pr due to their hybridization with the 4/^-shell of 
Pr. This view has been supported and extended by band 
structure calculations of Liechtenstein and Mazin, show- 
ing that with increasing radius of the RE just for |Efi-a| 
new band with 4/-character crosses the Fermi leveljH3'Ell 
and by Wang et al., who explicitly included a third band 
from the chain layers to describe the hole concentration, 
as a function of the Pr-concentration in the RE-layer.Eil 
As an alternative to the hybridization models it has been 
discussed that Pr might push theJaale states in a band 
associated with the apex oxygen.L3c3 

Experimentally, the hybridization shows (among other 
techniques) in detailed Pr-0 bond length analysisp3 in 
the stj-png broadening of the critstaJ field transitions in 
INSMta and x-ray absorption.LjO Using polarized x- 
rays Merz et al. were recently able to show that the ori- 
entation of the 2p-hole states in the structure is indeed 
different in Y123 and Prl2i giving very strong support 
to the hybridization modelEj Clearly, a microscopic char- 
acterization of the localized hole states is very desirable 
in order to distinguish between the models and address 
open problems such as the mechanism of the carrier lo- 
calization or of the role of the Pr- and Cu-magnetism at 
low temperature. 

In RE 123, doping of the Cu02-planes is achieved by 
variation of the oxygen content of the chain layer, and 
one of the basic ingredients of hybridization models is 
that this doping mechanism is independent of the RE- 
ion. A model for the dependence of the hole count on 
the oxygen concentration in the chains y has to address 
two problems, namely the distribution of oxygen in the 
chains, and ths hole count induced by a given oxygen 
configuration.Eil The non-random oxygen distribution in 
the chain layer is best described by the asymmetric next- 
nearest neighbor Ising modeLE^E^ It describes the various 
superstructures observed at intermediate y (most impor- 



tant the formation of chains in Y123 at y = 0.5) which 
depend strongly on the RE-ion. Liitgemeier et al. have 
shown that Cu(l)-(chain-)site ^iQE. is a suitable method 
to investigate this distributionJ^^'E^I but they did not in- 
clude Prl23 in their studies, so we analyse this case below 
in some detail. 

The dependence of the hole count on the oxygen con- 
figuration with the chain length as the most prominent 
parametenb^s been subject to detailed theoretical studies 
in Y123Ej'Ej On the experimental side the electric field 
gradient (EFG) at the Cu(l)-sites is a sensitive probe 
of the-charge balance between the planes and the chain 
layer.tZl Ohno et al. used this to estimate the number 
of mobile carriers iii the planes of metallic, fully oxi- 
dized Pr^Yi^^l23.El Below we use the Cu(l)-EFG at 
various oxygen contents to verify experimentally the as- 
sumption underlying the FR model, that the CuO-chains 
dope holes into the Cu02-RE-Cu02-layers, independent 
of the presence of Pr. 

One might expect that the special arrangement of 
charges in the hybridized state proposed by Fehrenbacher 
and Rice shows in the crystal field parameters of the 4/- 
shell of Pr. Unfortunately, the lines from crystal field 
transitions in INS are severly broadened, giving rise to 
large uncertainties in the symmetry as welLas in the ener- 
gies of the low lying crystal field levels.E3'E£l NMR is very 
sensitive especially to the symmetry of the. ground state 
and the position of the lowest excitations O In a previous 
work we reported the first measurements of the ^^^Py- 
resonance in Prl23 for y = lS3 Most remarkably, we 
found a different ground state symmetry for Pr"^"*" from 
the one used to describe the INS spectra. Moreover, the 
virtually nonmagnetic state of Pr at low temperatures 
which is without doubt what is observed in our NMR- 
measurements is in clear contrast to the antiferromag- 
netic Pr-order detected by ncutmn diffraction and the 
homogeneous suzeptibility.LJcd'ErE3 

Our interpretation of the NMR-dataJjia&JDeen strongly 
questioned because of these findings ,OciH3 but no con- 
sistent explanation for both, the NMR and the neutron 
data has been proposed up to now. Therefore we ex- 
tended the investigation of the Pr-resonance to crystals 
over the whole range y = ... 1 of oxygen concentrations, 
and to crystals with different concentrations in the Ba/Pr 
solid-solution system. On this basis we confirm our ini- 
tial assignment of the Pr-resonance to Pr on regular RE- 
sites. In view of the strong doubts our brief description 
of the 4/— ground state in refill did encounter we give 
below a rather detailed discussion of our determination 
of the crystal field ground states of Pr in this structure. 
This leads us to a consistent description of the NMR and 
neutron experiments, basically by assuming that the in- 
terpretation of both has been correct, and assigning two 
different electronic states of Pr to the conflicting experi- 
mental observations. 

Finally we discuss the impact of this model on the 
understanding of the low temperature magnetism of 
Pr and the other rare earths in this remarkable ma- 
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terial. As indicated above, the Pr-magnetism with its 
Neel-tempjpatui|e-| up to 18 K in Prl23 and in related 
cupratesc3E3cj e3 is as outstanding as the suppression of 
superconductivity. It is tempting to relate the two effects 
by the argument that the 4/ — 2p7r— hybridization sup- 
presses superconductivity by switching on pair breaking 
at Pr-moments, and simultaneously induces the high TV. 
via a substantial enhancement of the 4/ — 4/-exchange.El 
We would like to caution about this argument: While 
there is ample evidence for the hybridization itself and 
we present more below, the contribution of pair break- 
ing to the suppression of superconductivity is much less 
clear (see above), and the role of the 4/— exchange in 
Pr-magnetic order seems even more open. In fact the 
exchange is known to be small for all other rare earths, 
the ordering of the other 4/— moments is mostly due to 
dipolar interactions. l2I This is most convincingly demon- 
strated by the small dependence of on magnetic di- 
lutioiB-ituthe rare earth sublattice for all magnetic RE- 
ions.E3Til All magnetic RE, including the exceptional 
singlet ground-state Pr, order well above the 40% per- 
colation threshold for dilution of the square lattice with 
nearest neighbor exchange. The order is, therefore, as- 
cribed to the long range dipolar coupling of the magnetic 
rare earths, which is clearly not an option in the case 
of Prl23. The dependence of Tm on oxygen content is 
also puzzling, if hybridization enhanced exchange is the 
origin. At y = 0, where there is no hybridization of hole 
states, Tjv is still as high as 12 K. In Prl23, Tm increases, 
with the hole count, while it decreases for Nd in Ndl23.Ea 
In the concluding section of the discussion we show that 
the increase of Tjq of Pr is even more puzzling in view 
of our NMR-data, and discuss induced magnetism of the 
Pr-sublattice as an alternative model. 



II. EXPERIMENTAL DETAILS AND RESULTS 
A. Sample characterization 

The main diffculties in the experimental verification of 
the models arise from the complex doping mechanism in 
the RE123 series on the hand, and from the strong in- 
fluence of crystal preparation on the material properties 
on the other. A remarkable example is the preparation 
of crystals under special conditions that do become su- 
perconducting, with Tcmax even highM:-than in the other 
RE123 (> 100 K under pressure) H'E^ In order to pro- 
vide a reliable basis for comparison of results obtained 
with other samples we present in this subsection a rather 
detailed description of the sample preparation. 

High quality Prl23 crystals of a size suitable for NMR 
(> 10 mg) are still very difficult to prepare. One ma- 
jor problem is the contamination of the Cu(l)-sublattice 
with Al from AI2O3 crucibles (up to 30%), because-this 
has a strong influence on the oxygen distribution^ as 
well as on the magnetic structure of the Cu(2)-sublattice 



at low temperatureEa We avoid the Al-contamination by 
use of MgO-crucibles, which results in a much smaller 
contamination by Mg. (« 1%) because of the higher melt- 
ing point of MgOEj The best stoichiometric crystals 
available are grown in BaZrOs-crucibles, but they have 
been too small for NMR up to now. 

A second problem is the existence of a finite solid so- 
lution range for the EE/Ba-sublattices for the light rare 
earths (Nd, Pr, La).EZl"EJ The solid solution range is a 
four-dimensional volume in the phase diagram spanned 
by temperature T, partial oxygen pressure p{02), and the 
(2D-) triangle of cation concentrations, sketched for fixed 
T and p{02) in fig.|l|. The stable compositions form the 
surface of this volume in the phase diagram with the el- 
lipses in fig.|l| (top) indicating cross sections. Also shown 
are the liquidus surfaces and the conodal lines connect- 
ing them to the solid solution systems at two different 
p{02)- In the absence of other metals which may sub- 
stitute for Cu (e.g. Al from the crucibles) one may as- 
sume a constant Cu content. In this case the ellipses 
in fig.l degenerate to horizontal lines at the concentra- 
tions Ri+a;Ba2-a:Cu306+j,, and the solid solution range 
has only the three degrees of freedom T, p{02)^ and x. 

The accessible range of x depends on temperature, oxy- 
gen partial pressure, and on the rare earth radius. In 
general, the RE-content x increases with increasing oxy- 
gen partial pressure and rare earth radius, and with de- 
creasing temperature, at least in the vicinity below the 
peritectic temperature, We grew Prl23 crystals by the 
slow cooling method.Lj The Pr/Ba ratio of the growing 
crystals was set by the Cu/Ba ratio of the flux via the 
conodal line to the solidus surface. Temperature and 
oxygen pressure for variation of x were chosen as marked 
by the dots in the phase diagram to obtain crystals with 
high and low Pr-content along the conodal lines indicated 
in the flgure (see also tab.|). 
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TABLE L Preparation parameters and characterisation of 
the crystals. The flux composition is given in the first block, 
followed by the preparation atmosphere and the temperature 
program, with a slow decrease at the indicated rate in the 
given interval, and a quench or slow cooling to RT. The mor- 
phology and possible growth twins are indicated in the last 
block. 
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FIG. 1. Top: Schematic phase diagram of the Pr-Ba-Cu-O 
system at 940° C in 1 bar O2 (solid lines) and in 63 mbar O2 
(dotted lines). Stoichiometric phases are indicated by filled 
squares. A few conodal lines connecting the liquidus surface 
with the Pr/Ba solid solution system (Prl23ss) are shown 
with the starting compositions for the crystal growth marked 
by the dots. Bottom: Schematic phase diagram along the 
horizontal line corresponding to Pri+i,Ba2-iCu307_y in the 
upper part, again for 1 bar O2 (solid), and 63 mbar O2 (or 
300 mbar air, dotted). 

The two crystals prepared to obtain off-stoichiometric 
X were oxidized in 1 bar flowing oxygen in the temper- 
ature range 600-300°C to obtain y = 1. The crystals 
of the series with various oxygen content were prepared 
along the route denoted a; = in the table, the oxygen 
content was set in a similar tempering step under suitable 
oxygen pressure. The oxygen concentration y expected 
from the tempering conditions is in all cases in accord 
with the relative NMR-line intensities of the Cu(l)-0- 
coordinations in the chains (see below). The error in y 
estimated from the NMR-spectra is 5%. We expect that 
we may not reach the end members (especially y = 0) hy 
approximately that margin, but will nevertheless denote 
them for brevity with y = and y = 1 respectively. We 
did not succeed to prepare superconducting Prl23 along 
this route, all samples are semiconducting. 

At this point we would like to note that the local 



Pr/Ba-ratio of a given sample with fixed overall x can 
be influenced by an additional heat treatment at higher 
temperature as well, but the composition of the resulting 
crystal is inhomogeneous. For small deviations from the 
equilibrium surface a crystal may decompose into a mem- 
ber of the solid solution system with a different RE/Ba 
ratio, and either other stable phases like BaCu02, or 
a second member of the solid solution system. Larger 
deviations from equilibrium may lead to a spinodal de- 
composition where, for example, the Pr/Ba-ratio varies 
laterally and changes with time. The final state of these 
decomposition reactions depends on the initial RE/Ba- 
ratio and whether or not this cation ratio is located inside 
the four-dimensional volume of the solid solution system. 
In the case of a spinodal decomposition the regions of dif- 
ferent Pr/Ba ratios within a sample are very close to each 
other (10 - 100 nm), and any integral chemical analysis 
will reproduce the composition of the untreated sample. 

Unfortunately, the substitution of Ba by Pr (abbrevi- 
ated [Prjsa below) is hard to detect by standard scat- 
tering techniques because of the similar scattering cross 
sections of Pr and Ba. We have, therefore, up to now 
in our preparation little quantitative control of x. A 
unique way to determine the [Prjsa-concentration is 
the magnetic signal in polarized neutron scattering from 
these defects. Ed Crystals prepared along the route de- 
noted X = were analyzed by Markvard et al. using 
this technique, and a [PrJ^a-concentration near 5% was 
detected£3 We therefore expect x ~ 0.05 in the crystals 
we tempered to obtain various oxygen contents. 

Neutron diffraction showed that Pr on Ba-sites is not 
responsible for the extraordinary suppression of super- 
conductivity in Prl23. Kramer et al. report that 
[Pr]Ba and [Alfflsa have similar influence on the materi- 
als properties .0 The RE-ion binds oxygen on otherwise 
empty neighboring 0(5)-oxygen sites betsueen chains, dis- 
rupting the chain order in the process.E3 The defect is 
expected to reduce the hole count in the planes and pro- 
vides a strong scattering center in the chain layer, where 
it may—CQiLtribute to suppress ID-conductivity in the 
chains.EETEj However, the experiments indicate that the 
influence of the defects is limited to the chains, an inde- 
pendent mechanism not active in Ndl23 must be present 
in Prl23 to suppress the MIT and superconductivity in 
the Cu02-planes.Ej We emphasize, however, that [Pr]sa 
may well influence the low temperature magnetic struc- 
ture of the Cu(2)-sublal:li£;e by inducing magnetic mo- 
ments on Cu(l)-sites.Lj'0 NMR and neutron diffrac- 
tion consistently showed that Cu(l)-moments induced 
by (nonmagnetic) Al couple the planes ferromagnetically 
along the c-axis and can induce a rotation transition, 
of the Cu(2)-structure to avoid Cu(l)-frustrationE3'L£l 
A similar magnetic transition in Al-free, semicmiduct- 
ing Ndl23 may well be due to Nd on Ba-sites,c30 and 
Rosov et al. reported ordered Cu(l)-moments in Prl23 
at low temperature.L3 Magnetic-pu(l) is not detected in 
Cu(l)-NQR (see below and ref.L3) and is, therefore, most 
probably localized at defect sites in the structure. 
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While X is always positive for Nd, the accessible range 
in the Pr/Ba solid solution system may extendptp the 
Ba-rich side, as is the case for La/Ba in Lal23.E3EZl Sub- 
stitution of Pr (+3) by Ba (+2) should increase the hole 
count in the planes, and in fact it has been suggested that 
the suKjeijcoiiductivity in Prl23 crystals prepared by Zou 
et aLLaS'tZLpught be due to heavy Ba-doping of the Pr- 
sublattice.BIIj Note that disorder on the Pr-sublattice at 
smaller defect concentrations has been invoked as a vi- 
able mechanism to localize the hole states in the 4/ — 2p7r- 
hybridization bandEj 

The symmetry and the morphology of the oxidized 
crystals depend in a characteristic way on x that we 
also find in Lal23 and Ndl23. Ba-rich and stoichio- 
metric crystals exhibit an orthorhombic structure with 
twins, whereas the Pr-rich crystal kept its tetragonal 
structure even after additional oxidation treatments un- 
der high oxygen pressures. This result was expected be- 
cause the extra oxygen ions on the 0(5)-sites between 
the Cu(l)0-chains reduce the orthorhombicity. The mor- 
phology of the crystals changes with increasing RE/Ba- 
ratio from isometric blocks, sometimes with additional 
(lOl)-faces, formed by stoichiometric samples, to the for- 
mation of (100)7(001) growth twins and, finally, even to 
(001)/(111) growth twins at large x. All orthorhombic 
crystals investigated in this work are twinned with re- 
spect to the a- and 6-axes, but we find no indication for 
c-axis twins in the NMR-spectra. 
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FIG. 2. Inverse low temperature a, fe— plane mass suscep- 
tibilities showing the low temperature magnetic transition. 
The Cu(2)-sublattice orders antiferromagnetically near room 
temperature.i— iBeiep Tni (arrows) Pr carries a moment of 
0.3 . . . 0.9/iB.LJCjLJ The transition is very sensitive to sam- 
ple quality, Uma et al. observed a transition in-two steps in 
high-purity crystals grown in BaZr03-crucibles.L3 

In fig.|^ the low temperature susceptibility in a field 
of 1 T perpendicular to the c— axis clearly displays the 
cusp associated with the reorientation transition in the 
Cu(2)-sublattice and the appearance of the Pr-moments 
at various oxygen concentrations. The size and shape of 
the feature is very sensitive to sample qualitit^and even 
depends on magnetic history of the sample. l3 In clean 



samples the transition temperature drops from near 17 K 
at 2/ = 1 to 12 K at 2/ = 0, similar to what we find in 
our crystals. The transitions occure at 14.9 K and 12 K 
for the Ba-rich and the Pr-rich sample respectively (not 
shown) . 

B. The Cu(l)-0-Chains 

Fig.^ shows the zero field spectra, hgJJ the field sweep 
spectra in Prl23 at various oxygen concentrations x. The 
spin-echo spectra were taken at 4.2 K and at comparably 
long pulse distances (w 40^s) to suppress the contribu- 
tion of the Cu(2)-sites. The frequency range of 19 to 33 
MHz for the NQR-spectra of the Cu(l)-sites in the an- 
tiferromagnetic RE123 (RE^^ Pr) and the assignment of 
the lines to the different ojwaea coordinations of Cu(l) 
is by now well established.oEj The line positions and 
the relative intensities are similar to the ones observed in 
the other RE123, so we can safely use the same assigns- 
ment (Cu(l)4 = 21.8 MHz, Cu(l)3 = 23.3 MHz, and 
Cu(l)2 = 30.1 MHz for ^^Cu). The subscripts denote 
the number of oxygen in the Cu(l)-coordination shell, 
four for Cu(l) in a full chain, three for a chain end po- 
sition, and two for the apex oxygen in an empty chain 
position. 

The lines are significantly broader than in the other 
RE123 systems, where widths well below 100 kHz are 
observed in the homogeneous end members. The broad- 
ening may be due to structural defects inducing a dis- 
tribution in the electric field gradients (EFG), and to 
inhomogeneities in the magnetic structure of the Cu(2)- 
lattice inducing inhomogeneous internal magnetic fields. 
The two cases can be distinguished from spectra obtained 
in a large external field, where the frequency of the cen- 
tral transition is to first order independent of the EFG, 
so its linewidth is dominated by inhomogeneities of in- 
ternal magnetic fields, while the distance between the 
satellites and the central ransition is determined by the 
component of the EFG-tensor along the magnetic field, 
so the satellites are broadened by both, distributions in 
the EFG and in the internal magnetic field. In fact, both 
cases have been observed, Nehrke et al. found a signifi- 
cant narrowing of the Cu(l)-lines above T7V2 at ?/ = 0, a 
clear indication that magnetic disorder is present at low 
temperatures, while Grevin et al. confirmed our result 
that the magnetic transition at Ttq2 does not infiuence 
the linewidth at ?/ = 1, instead they find that it i s du| e to 
the formation of a charge density wave at 100 KtH-Ea 

The fluctuations at low temperature confirm, howevec, 
that the transition at Ti<!i in our crystals is magnetic.tll 
The insets to fig.^ display the systematic change of 
this behavior with increasing oxygen content. The in- 
sets show the ^^Cu(l)-signal intensity, corrected for the 
Curie-law of nuclear magnetization, versus temperature, 
measured across T7V2 (arrows). The spin echo is mea- 
sured at low temperatures with a repetition time shorter 
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than the relaxation time Ti for equihbration of the nu- 
clear spin system. Under such circumstances one expects 
that the amplitude increases when the relaxation time 
decreases, as is expected near a magnetic phase transi- 
tion. The stoichiometric crystal at y = (top) clearly 
shows this behavior. Ti decreases with increasing tem- 
perature and a pronounced minimum of Ti at Tn2 leads 
to a peak in signal intensity. With oxygen doping Ti be- 
comes shorter in the low temperature phase, leading to a 
higher signal than above and the relaxation induced 
peak at the transition is lost, presumably swamped by 
the large background of fluctuating fields in the ordered 
phase. At y = 1 finally, no significant effects of the tran- 
sition could be detected in this simple way. However, we 
found evidence for two contributions to the fluctuations 
in this temperature range for the fully oxidized crystal 
in full Ti -measurements. One component is almost in- 
dependent of temperature and may be quadrupolar or 
magnetic, the other shows a broad minimum at|-2jv2 and 
is ascribed, therefore, to magnetic fluctuations .Ed 
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FIG. 3. Zero-field NQR spectra of the Cu(l)-sites in Prl23. 
The assignment of the lines to the two- three- and fourfold 
oxygen coordination of Cu(l) corresponds directly to the one 
for the other RE123.The insets show the temperature cor- 
rected echo amplitude at constant excitation conditions ver- 
sus temperature to illustrate the systematic change in the 
infiuence of the fiuctuations at Tn2 (arrows, see text). 
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FIG. 4. Field-sweep spectra of the Cu(l)-sites in Prl23 at 
81 MHz and 4.2 K with the field along the (twinned) a- and 
6-axes. Together with similar spectra with the field along the 
c-axis (not shown) the positions of the quadrupole satellite 
transitions determine the electric field gradient tensor, their 
intensities the distribution of oxygen in the chains. 

The point symmetry of the Cu(l)-sites implies that 
the principal axes of the EFG-tensor for all oxygen co- 
ordinations are along the crystallografic axes. The full 
tensor is, thereby, determined by the assignment of the 
eigenvalues T4, Vy, T4 (in ascending order of absolut val- 
ues) to the axes, by the size of the largest eigenvalue 14, 
and by the asymmetry parameter yy =| Vx — Vy \ /Vz- 
We determined these values from fits to the field-sweep 
spectra with the field along the c— axis, and along the 
a—, 6— axes. The values are compared in table || 

TABLE II. Comparison of the largest eigenvalue T4(in 
10~^^y/m^), the asymmetry parameter r;, and the eigen- 
value of the EFG along the c— axis, for the three 
Cu(l)-coordinations in Y123 and Prl23. 
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with the ones observed in Y123. We cannot distinguish 
between a— and 6— axis in our twinned crystals, therefore 
we give only the EFG-component along c. The tetragonal 
symmetry of Cu(l)2 implies = 0, as is observed. 

In fig-H we compare the field-sweep spin echo Cu(l)- 
spectra of the three fully oxidized crystals, again with 
the external field applied along the a-, and 6-axes. The 
fourfold coordination Cu(l)4 characterized by the EFG- 
tensor ||I4z|| ~ W^wW or 77 w 0.9 is the only one which 
we could identify in our crystals, as is expected for 
filled chains. All Cu(l)-spectra are centered at the Cu- 
resonance of the metal, showing that the Cu-chain sites 
detected in NMR carry no moment in these crystals. We 
emphasize that this means only that there is no homo- 
geneous magnetization on the Cu(l)-sublattice. Cu(l)- 
moments localized at defect sites in the structure may 
well be undectable by NMR due to short relaxation times. 
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FIG. 5. Comparison of the field-sweep spectra of Cu(l) at 
81 MHz in Prl23 at y — 1 prepared in different oxygen atmo- 
spheres (see tab.^. The vertical lines at the center transitions 
mark the Cu-resonance from the wire of the pick-up coil, and 
at the '''^Cu-satellites the position in the stoichiometric sam- 
ple (bottom). The full lines are simulated spectra (see text 
for parameters). For the stoichiometric sample the decompo- 
sition into four subspectra (two isotops in a field along the 
local X- or i/-axis of the EFG) is also indicated. 



The spectrum of the Pr-rich sample (top) shows a 
larger splitting than the stoichiometric one (bottom) , an 
inhomogeneous broadening of the satellite lines, but no 
significant broadening of the central line. The high field 
shoulder at 7.4 T of the ^'^Cu-central transition can be 
assigned to the (111) growth twins mentioned above. As 
discussed above, the broadening of only the satellite lines 
shows directly the presence of disorder in the charge dis- 
tribution of the chain layer, clear evidence for the sub- 
stitution of divalent Ba by tetravalent Pr and the dis- 
order this defect site induces in the oxygen chains. On 
the other hand, the unchanged central line width proves 
the absence of internal magnetic fields at the Cu(l)- 
site, indicating that the magnetic structure of the Cu(2)- 
sublattice is not affected by the defects. 

The spectrum of the nonstoichiometric Ba-rich crystal 
(center) is, in contrast, best described by two subspec- 
tra. One with app. 30% relative intensity is the same as 
for the stoichiometric crystal, accordingly the environ- 
ment of 1/3 of the Cu(l)-sites is unchanged. A signifi- 
cant broadening of the central transitions together with 
a severe broadening of the satellites shows the presence 
of inhomogeneous magnetic fields together with charge 
disorder for the majority of the Cu(l)-sites. The inho- 
mogeneous width of the central transition corresponds to 
internal magnetic fields of w 0.2 T, in full agreement with 
the transferred hyperfine fields from the Cu(2)-momcnts 
at the Cu(l)-site (w 0.1 T) that are known from stud^ 
ies of the so-called AF-II structure in Al-doped Y123.E3 
The defects in this crystal do, therefore, introduce in- 
homogeneities in the magnetic structure of the Cu(2)- 
sublattice. 

The most plausible explanation for the fact that the 
size of the internal fields agrees well with the known 
transferred hyperfine field from the Cu(2)-moments is 
that the divalent Ba dopes a hole state into the dx2-y2 — 
2pcr-band of the Cu02-planes, similar to the well known 
case of Ca-doping on the RE-site in RE123. One might 
expect that this will form at low Ba-concentrations in 
the RE-sublattice a localised Zhang- Rice singlet, which 
in turn would lead to a Cu(2)-moment missing in the 
configuration of the neighbouring Cu(l)-site. The intrin- 
sic superconductivity of Prl23 crystals prepared under 
special conditions discussed recently is in view of these 
results most probably due to these holes becoming mobile 
at higher Ba-concentrations, just as in the other high-Tc- 
cuprates. 

C. Pr-Resonance 

We observed Pr-resonances very similar to the ones re- 
ported in ref.ref. ^ in all crystals except one specimen 
with the nominal concentration y = 0. In a second re- 
duced crystal we found the signal, but its intensity is very 
small. Fig.^shows the field sweep spectra at y ~ 0.2, fig.0 
shows the gyromagnetic lines for all samples. The signal 
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FIG. 6. Field-sweep spectra of Pr in slightly doped 
{y = 0.2) Prl23 with the field along the c-axis (top), and 
the a- and b-axis (bottom). At frequencies below ~ 40 MHz 
there is increasing overlap with the Cu(l)-spectra, especially 
Cu(l)2 at small oxygen concentrations. 




4 5 6 7 8 
Pi„H„ [T] ti„H„ [T] 

FIG. 7. Field dependence of the Pr-resonance frequency 
with the field applied along the c-axis (right), and along the 
a- and b-axes (left) for the different oxygen concentrations 
y. The full lines were calculated for the van Vleck paramag- 
netism of the Pr-4/-shell in the crystal field as described in 
the discussion. See Nehrke et al. for a detailed discussion of 
the line splitting in small in-plane fields. 

is clearly due to ^''ipr since this is the only ion in the 
structure which may show van Vleck paramagnetism and 



the corresponding anisotropic enhanced gyromagnetic ra- 
tios 7e// (the slopes at high fields, see below). 

7e// is almost indendent of y, only for y — we find a 
higher value than in doped crystals with the field along 
the c— axis, and a relatively small one with the field in 
plane. The fully oxidized, stoichiometric crystal is the 
only one where an orthorhombic distortion is detected 
by a splitting of the Pr-resonance line at low external in 
plane fields (see ref. |l^ for a detailed discussion of this ef- 
fect). This does not prove, however, that the symmetry 
is tetragonal on a local scale for y < 0.7 in our crys- 
tals. A splitting similar to the one reported for y = 1 
might be hidden, because the smaller Pr-signal at lower 
y overlaps with the very complicated Cu(l)-spectrum. 
The linewidth is seen to be roughly proportional to the 
frequency or the optimum field, indicating that the inho- 
mogeneous broadening is due to a distribution of the local 
gyromagnetic ratios with a width of app. ±10%. There 
is no indication of any unresolved-|quadrupole splitting 
like the one observed in Pr2Cu04j^ presumably because 
the EFG is even smaller at the nearly cubic Pr-site in 
the Prl23 structure. We point out that the Pr-linewidths 
turn out to be independent of oxygen content, in contrast 
to the Cu(l)-resonances. 
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FIG. 8. Pr-field-sweep spectra at 102 MHz, 1.3 K, in a field 
-Boljc. The amplitudes are scaled to allow comparison of the 
line shapes. The signal in Pr~123 is smaller by more than 
one order of magnitude. 

Comparison of the Pr-spectra in the three fully oxi- 
dized crystals in figjs] shows immediately that the effec- 
tive gyromagnetic ratios and the linewidths of the sto- 
ichiometric and the Pr-rich sample are identical. The 
intensities of the two signals have been scaled to the 
same maximum for this comparison. The result shows 
most clearly that the local electronic state of the Pr- 
ions contributing to the signal is the same in the stoi- 
chiometric and the Pr-rich crystal. As indicated above, 
there is no additional broadening in the Pr-spectrum 
for off-stoichiomctric oxygen concentrations, in contrast 
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to the obvious inhomogeneous broadening of the Cu(l)- 
resonance. This gives strong support to our earher as- 
signment of this signal to Pr on regular RE-sites and not 
to the [Prjsci-defect, because the regular site is shielded 
from the structural disorder in the chain layer by the 
Cu02-planes. 

In contrast to the similarities of the Pr-resonance on 
the Pr-rich side of the solid solution system we have 
scarcely been able to detect any Pr-resonance in the Ba- 
rich sample. Again this supports our assignment of the 
Pr-signal, since Ba substituting for Pr introduces disor- 
der in the RE-sublattice and it has to be expected that 
large changes in the effective gyromagnetic ratios of the 
neighbouring Pr-sites will effectively wipe out the Pr- 
signal from that region. 

III. DISCUSSION 
A. Hole doping 

The hole doping mechanism of the Cu02-planes with 
increasing oxygen content of the chains in the RE123- 
structure must be analyzed in two separate steps. First, 
the distribution of oxygen on the chain sites has to be 
determined, then the number of holes induced by the 
various Cu(l)-0-configurations has to be found. Summa- 
tion over the cluster probabilities times the corresponding 
hole counts should then yield the overall hale concentra- 
tion for a given oxygen concentration yE3 Cu(l)-NQR 
may, in principle, contribute to both questions. The res- 
onance frequencies determine the quadrupolc splittings 
and reflect the charge distribution, the relative intensi- 
ties of the lines correspond to the probabilities of the 
different oxygen coordinations of Cu(l), so they give in- 
formation on the oxygen distribution. 

The EFG at Cu(l)-sites depends sensitively on the 
charge distribution on the chains. It is hard to calcu- 
late the EFG in a correlated electron system from first 
principles, but for the chain layer local density calcula- 
tions show convincing results.H The main contributions 
arise from the charge distribution within the first Wigner- 
Seitz cell. This justifies the phenomenological ansatz to 
separate the EFG-tensor into a contribution from the in- 
completely filled Cu-(ia;2_y2 orbital and one from the p- 
orbitals of the neighbouring oxygens. Ohno et al. use for 
the Cu(l)-site a Sd-contribution of 117 MHz/hole and 
ti^pv'q = —67 MHz for the contribution from n^p holes 
on the four oxygen neighborsS, making the quadrupole 
splitting very sensitive to the charge distribution in the 
chain layer. Nevertheless, from table|l|the EFG of Cu(l) 
in Y123 and Prl23 is very similar for all oxygen coordi- 
nations. Vz in Prl23 is systematically smaller by « 5%, 
which might well be due to the differences in the lattice 
constants, and rj of Cu(l)4 is somewhat reduced in Prl23, 
which might be connected to the relatively large defect 
concentrations. Furthermore, we found no evidence that 
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FIG. 9. Comparison of the mean Cu(l)-0-chain length 
versus y determined from the relative intensitieS|— of the 
Cu(l)-lines in RE123 with RE=La, Tm, Gd, Nd, Y P), and 
Pr (this work). The full line sketches the expected behav- 
ior for a random distribution of chain oxygen, the dashed for 
repulsive O-O interaction V2, and the dotted for attractive 
V2. The bottom part shows the corresponding distributions 
of oxygen on the chains at half filling. This dependence of the 
oxygen distribution on the RE-radius may lead to a smaller 
doping of the Cu02-layer at intermediate y for larger RE, but 
the influence of the RE-radius on Tcmax = 89 ... 96 K is only 
small. 

the EFG at any of the Cu(l)-coordinations depends on 
the oxygen content, in contrast to the expectations in 
case of a localization of the carriers in the ID-chains. 
We conclude that the results strongly support the ansatz 
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that the substitution of Y by Pr has very httle influence 
on the carrier concentration in the chain layer. 

Cu(l)-lines from sites in CuO-clusters of different sizes 
are not resolved, so a direct determination of the cluster 
probabilities from the Cu(l) line intensities is not pos- 
sible. Following the procedure described by Heinmaa et 
al. we may estimate, however, the mean chain length 
(n) = 1 + 2/4/73 from the relative intesities Ia/Is of 
the three- and fourfold coordination. Note that there 
is no model for the oxygen distribution on the chain 
sites involved in this evaluation of (n). The small size 
of the crystals introduces a relatively large error in the 
orientation which has a strong influence on the field 
sweep spectra especially for Cu(l)3. To check for consis- 
tency we compare the intensities with the concentration 
y = 1^/2 -\- 1^. If we consider only spectra which could 
be fitted consistently in this sense we arrive at the de- 
pendence of the chain length on y compared in fig.^ with 
the corresponding results from Liitgemeier et al. 

Systematic errors for (n) are due to the overlap of 
the broad lines and to the assumption that no Cu(l)- 
contributions are missing in the experimental spectrum. 
It is quite possible that part of the Cu(l)-sites cannot be 
detected in NMR, either due to very large line broaden- 
ing from an inhomogeneous EFG-distribution, or due to 
fast relaxation induced by slowly fluctuating on-site mo- 
ments. As indicated above an ordered moment has been 
assigned to the chain sites in Prl23 by neutron diffrac- 
tion. We can exclude the presence of a small moment ho- 
mogeneously distributed over the whole Cu(l)-sublattice. 
At y > 0.4 we find no evidence for magnetic line broad- 
ening, and for y = the magnetic broadening of « 0.2 T 
together with the Cu(2) hyperfine coupling constant in 
RE123 (w UAT/^ib) sets an upper limit of 0.02^b/Cu 
for the crystals with low oxygen concentrations. How- 
ever, we cannot exclude the existance of localized mo- 
ments on Cu(l)-sites neighbouring defects. We empha- 
size again that such defect induced moments determine 
most probably the low temperature magnetism of the 
Cu(2)-sublattice in Al-doped samples,EJ and they might 
well be present on Cu(l)-neighbours of Pr on Ba sites. 

Despite this uncertainty in the integrated intensity of 
the Cu(l)-spectrum, which is common to some extent 
to all RE123-compounds with light RE, our experimen- 
tal chain lengths compare favourably with the expecta- 
tion for the chain length in the asymmetric next nearest 
neighbour Ising (ASYNNNI) model with vanishing V2, 
the interaction of two oxygen nearest neighbours with 
Cu(l) in between (see bottom of fig.||). The large, repul- 
sive interaction Vi drives the formation of chains with 
diverging length at zy = 1 in the ortho-I structure. The 
formation of the ortho-II superstructure of alternating 
empty and long chains in Y123 at y = 0.5 is due to an 
attractive V2 in the case of a small RE. If V2 is strongly 
repulsive the so-called herringbone superstructure is sta- 
ble at y = 0.5, and (n) will be one up to that point. This 
case is apparently nearly realized with the large La-ion 
on the RE-site. Nd with a radius in between is well de- 



scribed by a random distribution of the oxygen on chains, 
as is Pr. The origin of this influence of the RE- ion on the 
interaction potential of oxygen in the ASYNNNI-model 
is not clear at this time, but Prl23 fits smoothly into the 
trend. This finding underlines the conclusion above that 
Pr has no extraordinary influence on the electronic state 
of the Cu(l)0-clusters. 

The concentration of hole states in the Cu02-Pr-Cu02- 
trilayer should then be at all oxygen concentrations sim- 
ilar to the one in Y123, namely below the w 0.35/unit 
cell in (superconducting) Y123 at y = 1. Cluster and 
band structure calculations of the chemical potential as 
a function of the chain length in Y123 indicate that only 
clusters of a length above approximately 3|-Qxygen sites 
can induce hole states in the Cu02 planesHj Therefore, 
the reduced average chain length in Prl23 may lead to 
a somewhat smaller hole concentration in Prl23 than in 
Y123 at the same y, at least if this threshold behaviour 
occurs in Prl23 as well. Besides this small difference 
we arrive at the conclusion that on average at most ev- 
ery third Pr localizes a hole within the Cu02-Pr-Cu02- 
planes of the structure. 

B. The Pr 4/ ground state 

From the above arguments we are convinced that we 
can use the Pr-resonance to probe the charge distribu- 
tion in the Cu02-Pr-Cu02-trilayer by its effect on the 
CEF of the Pr-4/-shell. In cases where the NMR-signal 
of non-Kramers ions like trivalent Tm or Pr can be ob- 
served, rather detailed information especially on thp low 
energy part of the CEF-Hamiltonian is obtained^. An 
experimental determination of the crystal field at Pr from 
NMR is the more desirable since the spectra obtained in 
Prl23 by INS are heavily broadened and the difhculties in 
their analysis lead to-a considerable scatter in the CEF- 
energy level schemesll3. A more extensive description of 
our analysis of the crystal field Hamiltonian Hcf than 
was possible in ref. |l^ seems to be in order in view of the 
controversal discussion initiated by the results obtained 
for the crystal with the highest oxygen concentration. 

As stated above, the anisotropic enhancement of the 
gyromagnetic ratio {^'^^'y = 13 MHz/T without van Vleck 
contributions) unambiguously identifies the signal as due 
to ^"^^Pr. It is, however, more difficult to identify the 
position of the Pr in the lattice, and it has been argued 
that our NMR-experiments might detect the signal from 
the [Prjsa-defect sites. 

There is strong experimental evidence against the 
assignment to the defects. As indicated above, the 
linewidth of the Pr-resonance is independent of any 
oxygen- or Ba-site disorder in the chains, though this is 
clearly reflected in the Cu(l)-spectra at intermediate y or 
off-stoichiometric Pr/Ba cation ratio. Second, figs. ^ and 
[lO| show in accord with the detailed crystal field analy- 
sis below that the van Vleck susceptibility tensor is very 
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nearly tetragonal, and almost independent of y. This 
shows a tetragonal symmetry of the crystal field (CEF), 
independent of the oxygen contents of the chains. A 
symmetry close to tetragonal is in full accord with the 
expectation for the regular Pr-site at all y, because the 
site is far from the chain layer and shielded by the Cu02- 
planes. The CEF of Pr on Ba-sites, on the other hand, is 
expected to reflect the orthorhombic symmetry at y = 1 , 
where the neighbouring chains are filled, but also at small 
y, because Pr"^"*" is known to bind oxygen on the neigh- 
bouring 0(5)-sites. In addition, AUenspach et al. showed 
that the CEF of [Ki^Ba is completely different from the 
one on the RE-siteJ^ while we show below that the CEF 
we find from our results and the one determined from 
INS are rather similar. The strong asymmetry of this 
configuration has been observed in the EFG of La on the 
Ba-site£3 
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FIG. 10. Dependence of the efFective gyromagnetic ratios 
of Pr (top), and the Pr signal amplitude normalized to sample 
volume (bottom) on the oxygen concentration. The error bars 
in the right figure indicate along y the uncertainty from the 
relative Cu(l)-line intensities, along the intensity axis they 
estimate the uncertainty due to the different filling factors of 
the NMR-coils. The sample denoted throughout the text by 
y = for convenience is almost certainly not fully reduced. 
In a second crystal with the same Cu(l)-spectra we observed 
no Pr-resonance at all. 

The connection between the anisotropic 7e/ / observed 
in fig.0 and the CEF can be understood in terms of a 
single-ion description of the 4/-shell. In ref. |l8| we derive 
the nuclear spin hamiltonian 
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J 



.a 



A 



gjfj-B 



A 



J 



{M^)spla (3.1) 



for the nuclear spin I coupling with a gyromagnetic con- 
stant 7/27r = 13.0 MHz/T to the external field Hq, and 
via the hyperfine coupling constant Aj/h = 1093 MHz to 
the moment (M) in a 4f-shcll with angular momentum 
J and Landc- factor gj = 0.8. Xa ~ dMa/dHa is the 
static van Vleck susceptibility with the field along the 
crystallografic axis a = a,b, c, {Ma) sp is the spontaneous 
Pr-moment along one axis. The diagonalization of 71 for 
H[)\\M.sp or in the case of large external fields is trivial 
in the notation with the effective gyromagnetic ratio 7^ 
in the second part of eq. 3.1. The gyromagnetic ratio 7 
changes to 



7a = 7 + 



A, 



(3.2) 



which is independent of temperature and field as long 
as this is the case for the van Vleck susceptibility. This 
holds for temperatures and fields ksT, /io-ffo(M)ct ^ A^, 
where Aq, is the smallest CEF-splitting of 4/-levels con- 
nected by Ja- From the linear field dependence of the 
resonance frequency in fig.|^ it is clear that 7^ at 1.3 K 
is indeed independent of the field up to at least 9 T (5 
T for the field in the plane). The offset and the slope of 
the gyromagnetic lines determine immediately the spon- 
taneous Pr-moment (M)sp and the local susceptibility 
Xa respectively. 

The frequency offset is clearly smaller than 25 MHz 
for all samples. It corresponds to a very small or- 
dered moment of below hv/Aj — 0.023/iB/Pr in all 
cases, more than one order of magnitude smaller than 
the 0.3...1.2/XB detected with neutron diffraction or 
Mossbauer spectroscopy. It was perhaps this property 
which gave rise to the most severe doubts about the in- 
terpretation of our spectra of the fully oxidized sample. 

In order to understand the origin of this discrepancy it 
is usefuU to discuss first the CEF. We compare the three 
experimental components of Xa 



Xo 



A, 



(7a - 7) 



(3.3) 



with the ones calculated in the eigenbasis of a suitably 
modelled crystal field Hamiltonian Ti^f {Ti-cf \ k) = 
Sk I k)): 

Xa — 



keT 

k " I 
Y.\{k\Ja\l)\ 



' ]. (3.4) 

£1 - £k 



Because of the exponential factors in eq.3.4 the suscep- 



tibility is dominated by the matrix element of the total 
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angular momentum with the ground state and the low- 
est excited state with nonvanishing matrix element. The 
most convenient representation of Hcf for a calculation 
of Xa is by means of the Stevens equivalent operators O™ 
with crystal field parameters bim'- 



(3.5) 



It describes the two 4/-electrons in Russcl-Saunde^ 
coupling. The LiS-admixture from spin-orbit coupling,E3 
which is 3100 KJ23 and the interconfiguration coupling, 
i.e. the admixture of states fpspi the multiplets J = 5, 6 
considered by Hilscher et al.,t^ are neglected. This may 
be justified as a first approximation since the spin-orbit 
coupling is significantly larger than the crystal field split- 
ting. Then J = L -|- S is a good quantum number and 
the ground state of an isolated Pr^+ is the ninefold ^H4- 
multiplet, which is further split by Hcf, or by a magnetic 
field. The representation is based on the Wigner-Eckhard 
theorem, which guarantees that the matrix elements of 
Hcf, taken within a subspace of states | LSJMj) with 
fixed quantum numbers L, S, and J, are proportional to 
a suitable polynomial 0™ of angular momentum opera- 
tors. The equivalent operators and the numerical factors 
occurring in ihe Wigner-Eckard theorem are tabulated 
by HutchingsS 

The advantage of the equivalent operators 0[" is that 
the matrix elements can be calculated in a physically 
transparent way, since they contain only angular momen- 
tum operators, and usually a good fit of the energy spec- 
trum is possible. We follow, therefore, the common prac- 
tice in NMR and use below equivalent operators. Note 
that A J and gj enter both experimental parameters, the 
slope and the offset of the gyromagnetic lines. The val- 
ues for Aj and gj introduced above are the ones for the 
•^ff4-multiplet, and they allow for a fit to the gyromag- 
netic lines with reasonable parameters. This is not pos- 
sible if Aj is chosen significantly higher (to allow for a 
higher M^p), because then unphysically large moments 
result in high field. 

A good fit of the CEF energy spectrum from within 
the J = 4-multiplet alone does, however, not necessarily 
imply the same precision for the crystal field parame- 
ters. The analysis of the INS-spectra is usually based 
on a irreducible representation of Hcf by spherical har- 
monics and shows admixture of states from higher con- 
figurations. The conversion factors for the crystal field 
parameters Bi„i used in the irreducible representations 
and the him above may be determined by comparison of 
the coefficients witli_equal symmetry and have been tab- 
ulated by KassmarO. A comparison of Hcf in the two 
representations is difficult if the interconfiguration mix- 
ing has significant influence. In this case differences may 
occure in the eigenvalues of Hcf calculated from the irre- 
ducible representation and the ones calculated from the 
equivalent operator representation after conversion of the 
parameters. 



For 4/-electrons spherical harmonics up to the order 
Z = 2 • 3 = 6 and m < I can contribute to the series 
in eq.3.5. The point symmetry of Pr allows only for 
even l{— 0,2,4,6) and ■m{— 0,2,4,6 < I) in orthorhom- 
bic point symmetry, and the fourfold c— axis in tetrag- 
onal symmetry reduces this further to / — 2,4,6 and 
m = 0,A < I. Obviously it is impossible to determine the 
full set from the three components of the local susceptibil- 
ity Xa, and in Prl23 even the more detailed information 
from INS is insufficient. Fortunately, both methods are 
complementary: NMR is sensitive to the ground state 
symmetry and the splitting of only the low lying lev- 
els, while the INS-spectra suffer especially at low energy 
from the severe broadening, but well resolved lines are 
observed w 500 K above the ground state. 

A simultaneous fit of Hcf to the results of both meth- 
ods is viable only for the samples with high oxygen con- 
centrations y. The difficulty is that our ground state sym- 
metry for all y must be tetragonal with at most a very 
small orthorhombic distortion. The severe line broaden- 
ing of the INS-spectra does not allow to determine the 
ground state symmetry unambiguously for y = 1, but for 
y = any description with tetragonal symmetry is very 
poor. We have, therefore, apparently conflicting results 
from the two methods at y — 0. 

In order to solve the problem we propose to take the 
results of both experiments at face value and assume 
the presence of Pr in two different symmetries on the 
RE-sites. We emphasize that the presence of two dis- 
tinguishable electronic conflgurations of Pr was already 
proposed by Fehrenbacher and Rice, so we use below their 
nomenclature for the two sites. Here we only make use 
of the fact that the CEF should be different for Pr^+, 
the site without a hole state in its 4/-shell and in the 
2p7r-orbitals of the eight oxygen ligands, and for Pr^^, 
the site where a hole is localized in the 4/-shell and 
the 2p7r-ligand orbitals. The concentration of the Pr^^- 
conflguration is according to Fehrenbacher and Rice di- 
rectly the hole concentration rih determined above from 
the Cu(l)-resonance, that is at most « 0.35 holes per 
unit cell or roughly every third Pr-site. 

First we assign the Pr-resonance detected in NMR to 
the Fehrenbacher-Rice state Pr^^, and Pr'^^ to the signal 
obtained in INS (and most other techniques). From the 
NMR point of view this is natural because the intensity of 
Pr-resonance increases with y, that is, with the number 
of holes (see figjl^) . From the INS-results the assignment 
is necessary because a signal is present at ?/ = and the 
CEF observed is at all y similar to the one determined 
with very high accuracy for other RE in RE123, where 
the hybridization state RE^^ does not exist. 

The contradicting results of INS and NMR for the re- 
duced samples {y = 0) are now easy to understand. NMR 
observes a small signal from a few percent of Pr in the 
state Pr^^, which may be induced e.g. by a small amount 
of oxygen remaining on chain sites. In INS the small 
signal corresponding to these sites is probably not de- 
tected. On the other hand, the Pr-resonance from the 
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state Pr^+ dominating the INS-spectra may well be not 
observable in NMR, either due to fast relaxation, or due 
to a large inhomogeneous broadening from the EFG of 
the distorted 4/-shell and the magnetic hyperfine fields 
induced by 4/-moments. 

Near y = 1, where app. 1/3 of Pr is expected to be 
in the state Pr^^, they certainly contribute to the INS- 
spectra. Then we must seek for a simultaneous fit of Ticf 
to the NMR-data and the INS-spectra, weighted for their 
high sensitivity at low and high energies in the spectrum 
of TLcf respectively. Any differences in the CEF-energies 
of the two Pr-configurations must be compatible with the 
supposedly inhomogeneous linewidth of the INS-spectra 
at energies above app. 500 K. 

The method we use to fit of Ticf to the NMR-data and 
the high energy part of the INS-spectrum is motivated by 
the almost cubic oxygen coordination of the RE-site in 
RE123. We start from a CEF with cubic symmetry. We 
then improve this solution with a tetragonal distortion, 
which we may 'polish up' with an additional orthorhom- 
bic distortion and a local magnetic field. 
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FIG. 11. Determination of the crystal field parameters bim 
by successively lowering the symmetry in the fit from cubic 
(left), to tetragonal (center), to orthorhombic (right). The 
horizontal lines represent the crystal field splittings deter- 
mined by INS (left and top figures) and the van Vleck mo- 
ments of Pr determined from our NMR-data (bottom figures) . 
Note that the upper right figure contains only the levels at 
high energies, showing that a reasonable fit to the energies 
in the resolved part of the INS-spectra is achieved with this 
simple procedure. 

The three steps are shown in fig.|ll|. There is a general 
consensus from all experimental techniques that the cubic 
F^-triplet is the subspace underlying the quasitriplet at 
low energies of Hcf in Prl23. We therefore adjust the 
two parameters 640 , 650 of the cubic Hcf (see ref. ^ to 
fit the splitting between the ground state and the other 
levels, and the one between the other levels roughly to 
the experimental energies from INS. The left panel shows 
the regime of bgo for cubic symmetry, where the Fs-triplet 
is the ground state (640 is fixed at 0.335 K). We choose 
^60 = 0.001 K as a cubic approximation. 

Next we consider a tetragonal distortion by variation 



of beo with 654 = —0.021 K fixed at its cubic value (center 
panels). The F^-triplet splits into a doublet and a singlet, 
which are the ground state for bgo larger or smaller than 
the cubic value respectively (top part). The local sus- 
ceptibility (bottom part) Xz fits the local susceptibility 
observed in the NMR-experiment (horizontal lines) only 
if the singlet is the ground state. Xz is too large by two 
orders of magnitude if fego is larger than its cubic value 
(in the right half of the central panels). With the sin- 
glet ground state the tetragonal Xx (o) falls between the 
orthorhombic in plane susceptibilities at 650 = —0.002 
K. Simultaneously the splitting of the upper levels is de- 
scribed more or less correctly (top). At small y, where 
we find tetragonal symmetry, this procedure leads to the 
parameters given in tab. Ill, 

The crystal y = 1 is the only one where we found 
an orthorhombic splitting, and here we introduce the or- 
thorhombic distortion 622 shown in the right panel of 
figjl^. The changes in the high energy eigenvalues are 
small (top, note the suppressed zero), and even the eigen- 
vectors of Hcf change only by a few percent. The fit 
to the transitions observed in INS (horizontal lines) is 
already reasonable. The local susceptibilities (bottom) 
are well described at 622 = 1.9 K. We may now use the 
other parameters to 'polish the fit up', but we emphasize 
that this is far from unambiguous and any improvement 
is probably beyond the accuracy of the description with 
Stevens equivalent operators. The parameters given for 
y — 1 in tab. [II are an example for an almost exact fit to 
the local susceptibility and the energies above 500 K. 

TABLE III. Top: Crystal field parameters determined 
from the simultaneous fit to the high energy INS-spectra and 
to the local susceptibilities measured by NMR in the case of 
Prl23, X — 1. Since we find tetragonal symmetry for the local 
susceptibility we fit only to Xa for Prl23, y = 0, and searched 
for the nearest solution in the parameter space to the one of 
y = 1. Bottomirfligenvalues of TCcf as determined by Hilscher 
et al. from INStJ, and from our NMR-data. Values in brack- 
ets are not directly observed. Note that we did not intend to 
fit the energies in the case y — 0. 
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17 38 731 786 (825) (854) 963 1130 
132 142 726 (775) (832) (922) (956) (1126) 


Ecf(x^O) 
INS [K] 
NMR [K] 
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FIG. 12. Polar plot of the shape of the 4/^-wavefunction 
of Pr as determined from our NMR-data (top), and from 
INS (bottom). The (tetragonal) c-axis points to the top, the 
eight lobes along the local (lll)-directions point to the oxy- 
gen neighbours. Only one of the two possible orthorhombic 
ground states for INS is shown, the other one is rotated by 
7r/2 around the c-axis. 

Fig.|l^ shows a polar diagram of the charge distribu- 
tion in the two-electron ground states calculated from 
our crystal field parameters for y = 1, and from the con- 
verted crystal field parameters of Hilscher et al.li3 The 
state Pr^^ observed in NMR is tetragonal, the small 
orthorhombic distortion is not visible. The eight lobes 
pointing along the space diagonals of (nearly) a cube with 
oxygen at the corners resemble cicely the 4/-state en- 
visaged by Fehrenbacher and RiceE3. Note that there are 
four additional lobes in the plane perpendicular to the 
c-axis pointing to the edges of the cube. The ground 



state Pr'^+ observed in INS looks very different at first 
sight, but from the above discussion the connection be- 
tween both states is easy to visualize: Rotations of Pr^^ 
by 90 deg around the local x- or j/-axis (_L c) lead to the 
two remaining states of the cubic F'^-triplet, the doublet 
in the tetragonal symmetry. In a suitable orientation of 
the local x- and y-axes such a rotated Pr-^^-state resem- 
bles closely the state Pr'^+ in the figure, only the four 
lobes which were formerly perpendicular to c are more 
pronounced. Note that the orientation of the local x- 
and y-axes or these four dominant lobes is determined 
only when e.g. a Jahn- Teller distortion lifts the twofold 
degeneracy of this ground state in tetragonal symmetry. 

The physical origin for the different ground state sym- 
metries is in our model the one missing electron charge 
in the 2p7r-ligand orbitals of the oxygen coordination of 
Pr in the state Pr^^. In the presence of this hole the 
tetragonal CEF is on the left hand side in fig.|ll| with the 
singlet ground state. If the orbitals are occupied the CEF 
changes to the right hand side and the doublet becomes 
the ground state. The doublet degeneracy may be lifted 
by a Jahn- Teller distortion leading to the ground state 
observed by INS. 



C. The low temperature magnetic structure 

We now turn to the low temperature magnetic state of 
Pr. All crystals investigated here show the well known 
kink in the low temperature susceptibility, indicating a 
magnetic transition at temperatures increasing with oxy- 
gen content from « lOK to 18K (see fig.||). In the dis- 
cussion following our first observation of nonmagnetic Pr 
it became clear that there is a simultaneous transition 
in the magnetic structure of the Cu(2)-sublattice, but 
even taking this complication into account moments of 
at least O.S/is/Pr are required to describe the magnetic 
Bragg peaks at all y. 

The observation of a Pr-moment is consistent with the 
interpretation presented above if the state Pr'^'*' carries 
this moment, and the existence of two magnetically dif- 
ferent sites has been considered much earlier as a possL 
ble explanation for the broad ^""^Pr-Mossbauer spectra.EHl 
The fact that T/v increases with the concentration of the 
nonmagnetic site Pr^^ is, however, in clear contradiction 
to the concept that Pr-Pr exchange is the origin of the 
transition. Even the fact that the transition occures at all 
would be surprising at y = 1, because the magnetism of 
the RE-sublattice is known to be strongly 2-dimensional, 
and the percolation threshold for the square lattice with 
nearest neighbor exchange is only 0.5928. The 30% non- 
magnetic impurities introduced by hole doping would 
bring the lattice almost to the percolation. Obviously a 
similar difficulty with the exchange coupling arises with 
the small dependence of the Pr- T^r on magnetic dilution 
by other RE-ions. 

We propose, therefore, that the Pr-moments arc in- 
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duced by the transition in the Cu(2)-sublattice. The 
mechaaism is the same we introduced in our previous 
work,ll3 but now we have to consider two different Pr- 
sites. The different magnetic properties of the two Pr- 
sites are a direct consequence of the two ground state 
symmetries discussed above: The tetragonal symmetry 
of the state Pr^^ has as a consequence the cancehation of 
all transferred hyperfine fields from Cu(2) at the Pr-site, 
therefore the Pr-moment vanishes. The small moment 
indicated by the offset in the resonance frequency may 
be induced by dipolar fields from the Cu(2)-moments, 
which do not fully cancel by symmetry. 

On the other hand, the four pronounced lobes of the 
state Pr'^+ can point to the centre of the edges of the cube 
spanned by the oxygen, that is to four of the Cu(2)-sites 
(see flgjl^). Note that these moments are parallel in the 
structure proposed as a first approximation in ref. |l^ and 
sketched in the figure, and that they do have a parallel 
component in the refined structure determined later by 
neutron diffraction. In such a case the transferred hy- 
perfine fields at the Pr-site do not cancel. The size of 
the induced moments depends exponentially on the small 
splitting of the quasitriplet in the CEF (see eq.3.4). De- 
spite the considerable scatter in the CEF-parameters de- 
termined for Pr'^+ from INS it is clear that the splitting 
is smaller by at least a factor of two than the 130 K we 
obtained above for the state Pr^^, with the result that a 
transferred hyperfine field in the order of one Tesla can 
be sufficient to induce the observed moments. 

Note that no extraordinarily high Pr-Pr exchange is re- 
quired to explain the surprisingly high T/v in this model, 
therefore the problem of the small dependence of the 
transition temperature on magnetic dilution does not oc- 
cure. However, an exchange interaction of the order of 
1 K may help to understand at least qualitatively the 
peculiar orientation of the Pr-moments, halfway between 
the Cu02-plane and the c-axis.EiHi2l If the Pr-moments 
are induced by internal fields rather than due to 4/ — 4/ 
exchange interactions the magnetic structure of the Pr- 
sublattice is expected to reflect the one of the Cu(2)- 
sublattice, a feature which might be used to check our 
model. In this context it is interesting to note that Pr 
orders magnetically at exceptionally high Tjv in a number 
of related cuprates when it is situated, between antiferro- 
magnetic (insulating) Cu02-sheets.c3cZl To our knowl- 
edge there has been only one ijeeort of coexistence of 
Pr-order and superconductivityj23 but the sample was 
inhomogeneous and only partially superconducting, so 
magnetic order and superconductivity might be located 
in different volume fractions. 



with the one in Y123 confirms the assumption of the 
hybridization model, that the hole states doped with in- 
creasing y into the structure are not localized in the chain 
layer but in the Cu02-Pr-Cu02-trilayer. The relative in- 
tensities of the resonances corresponding to the different 
oxygen coordinations of Cu(l) in the chains show that 
there is no tendency to form 0-Cu(l)-0-chains with Pr 
on the RE-site, a result which fits nicely into the trend 
of the interaction of oxygen on chain sites with the RE- 
radius. 

Our investigations of the Cu(l)- and Pr-rcsonance in 
off-stoichiometric crystals in the Pr/Ba solid solution sys- 
tem show that the accessible concentration range for x 
in Pri_|_a;Ba2-2:Cu306+j, can be extended to negative x 
by preparation at reduced partial oxygen pressure. The 
fingerprint of the [_Ba]pr-defect in the Ba-rich crystal is 
a wipe-out of the Pr-signal and an inhomogeneous inter- 
nal magnetic field at the Cu(l)-site, which we ascribe to 
localized singlet states in the antifcrromagnetic Cu02- 
layers. [Pr\Ba in the Pr-rich crystal, on the other hand, 
does not infiuence the Pr-resonance or induce disorder 
in the magnetic Cu(2)-structure, but the structural and 
charge disorder in the chain layer is evident from the in- 
homogeneity of the Cu(l)-quadrupole splittings. These 
findings give strong support to our earlier assignment of 
the Pr-signal to regular (RE-) sites. 

The Pr-resonance was observed in all crystals except 
one fully reduced, stoichiometric specimen. In order to 
give a consistent interpretation of this NMR-signal as 
well as for the observations from INS and the volume 
probes, we propose that the two electronic Pr configura- 
tions introduced by Fehrenbacher and Rice correspond to 
two different ground state symmetries of Pr in the RE- 
layer. We provide evidence that NMR detects only the 
state with a hole localized in its 4/-shell and the 2p7r- or- 
bitals of the surrounding oxygen ligands, while the prop- 
erties observed with space-integral techniques like INS 
are dominated by the response from sites without such a 
localized hole. The model is supported by our count of 
localized carriers with the Cu(l)-resonance, the decreas- 
ing intensity of the Pr-resonance with decreasing oxygen 
concentration y, by our detailed analysis of the CEF of 
Pr, and it explains the discrepancy in the magnetic states 
of Pr observed in NMR on one hand, and INS and sus- 
ceptibility on the other. Finally, we argue against the 
assumption of a strong, enhanced 4/ — 4/ exchange as 
the origin of the extraordinarily high Pr ordering temper- 
ature and the suppression of superconductivity. As an 
alternative we propose that the Pr-moments are induced 
by internal fields which appear due to the reorientation 
transition of the Cu(2)-moments. 



IV. CONCLUSIONS 



We reported NMR/NQR investigations of the charge 
distribution in Pri_|_a;Ba2-xCu306+j, crystals at low tem- 
perature. Comparison of the EEC at Cu(l) in the chains 
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